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Summary 


Aerodynamic equations with nonlinear unsteady effects were formulated for an aircraft in a one- 
degree-of-freedom large amplitude motion about each of its body axes. The corresponding 
aerodynamic models were expressed in the form of indicial functions. The model formulation 
separated the resulting aerodynamic forces and moments into static terms, purely-rotary terms and 
unsteady terms. For model identification from experimental data it was assumed that the static and 
purely-rotary terms were known. The model identification procedure developed combines a stepwise 
regression and maximum likelihood estimation in a two-stage optimization algorithm which can 
identify the unsteady term and also rotary term if necessary. 

The identification scheme was applied to oscillatory data in pitch for two examples. The first 
example used the simulated data of a tailless aircraft, the second wind tunnel oscillatory data of the F- 
16XL aircraft. The results from both examples indicated that the two-stage optimization algorithm 
can converge to maximum likelihood estimates. The identified model from experimental data fit the 
data well, however, the accuracy of some of the estimated parameters was rather low, around 10%. 
The identified model was a good predictor for oscillatory data and data with ramp input. 


Symbols 

Aj,BjJ = 1,2,... 
a(a), b\ (a) 

a j'J ~ 0,1,2,... 
b 

h 

b[.,j = 0,1,2,... 

C a 

C L , C m 


spline terms 
polynomials in a 
coefficients in a(a) 
wing span, m 

indicial function parameter, sec -1 
coefficients in £q(a), sec -1 
general aerodynamic coefficient 
lift and pitching-moment coefficient 


1 



indicial functions 


C a ct (0’ C a q (O’ 

Ca, p (0’C ap (0> 

C Or (0 

c mean aerodynamic chord, m 

cj,j = 0,l,...,m coefficients in h{t\a ) 


D 

F(t; a) 

F a a (0’\(0’ 

F a q (0’%(0 

h(t;a) 

k 

£ 

m 

n 


P, q,r 


t 

V 

v(0 

y 

z 

a 

P 

A 

£ 

i 


operator, D = d/dt, sec 1 
model of deficiency function 

deficiency functions 

sum of exponential functions, eq. (5) 

reduced frequency, k = wl/V 

characteristic length, t = c/2 or i = b/l , m 

number of exponentials in h(t;a) 

number of data points 

roll, pitch and yaw rate, rad/sec 

time, sec 

airspeed, m/sec 

measurement noise at time (i - 1) A t 
variable defined by eq. (12) 
variable defined by eq. (25) 
angle of attack, rad or deg 
sideslip angle, rad 
increment 
equation error 
dummy integration variable 
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X 

<t> 

(0 


time delay, sec 

roll angle, rad 

yaw angle, rad 

angular frequency, rad/sec 


Superscript: 

time derivative 


Subscript: 

E measured value 


o 


nominal value 


Abbreviation: 


ML 


maximum likelihood 


SR stepwise regression 

SNR signal to noise ratio 


Aerodynamic derivatives: 

C A a ( cc ) = C A a ='^= L ’ for A = D,L,YJ, m , or n 

a “da 


for a : = p,q,r,a,a,/3 or /? 


_ pb qc rb ac n 6b 

and for a = • > a > ’P’ or 

auu ioi 2V 2V2V 2V 2V 


Introduction 

One of the first attempts to obtain unsteady aerodynamic characteristics of an aircraft from 
experimental data was reported in reference 1. Aerodynamic models included additional state 
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variables to those defining aircraft motion. These additional variables were used in defining unsteady 
effects. Experimental data came from wind tunnel and flight tests, parameter estimation used the 
ordinary least squares. Further improvements to modeling and parameter estimation procedures 
followed and are presented in references 2 and 3. Similar approaches to model formulation and data 
analysis by other authors can be found in references 4 and 5. 

References 6 to 8 present formulations of aerodynamic model equations in terms of indicial 
functions. The first of these references includes a method based on Fourier analysis of wind tunnel 
data from large amplitude oscillatory motion and motion generated by a ramp input. Estimation of 
parameters in references 7 and 8 was limited to models with linear aerodynamics. A different 
approach from the previous two is given in reference 9. The aerodynamic coefficients are specified 
as nonlinear functions of the motion variable and its rate of change. At the same time all the 
parameters in the model are considered as functions of frequency. 

In this report, the approach of references 7 and 8 towards modeling and parameter estimation is 
extended to cases with nonlinear unsteady aerodynamics. After the introduction, the report presents a 
development of mathematical models of an aircraft performing a one degree-of-freedom motion 
about one of the body axes. The models developed are then used in parameter estimation with 
simulated and real wind tunnel data from oscillatory tests in pitch. The problem of selection or 
determination of a specific model structure prior to parameter estimation is also discussed. The 
estimation methods are based on the least squares and maximum likelihood principles. Final models 
are assessed as to their ability to fit the measured data and predict the aircraft motion. The report is 
completed by concluding remarks. 

Postulated Models 

In this section mathematical models of an aircraft performing a one degree-of-freedom (d.o.f.) 
motion about each of the three body axes will be developed. These models will be applicable to 
aircraft harmonic motion, response to a ramp input or any other form of a single input. 
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Motion in Pitch 


For a one d.o.f. motion in pitch, the fundamental relations for the drag, lift and pitching moment 
are 


c a ( r ) - C a (a(r ),q(t )), a = D, L, or m 

Using the results of reference 10, each of the aerodynamic coefficients can be formulated as 


( 1 ) 


QM = Q(°) + J C a a (* “ U «( T )><?( T )) «( T ) dr 
0 
t 

+ yf Ca J l ~ T; a ( T )M T )) ?( T ) dT 

0 

where C a ( 0) is the value of the coefficient at initial steady-state conditions, C Ua (t -x; a(x),<y(x)) and 


C a (f-x; a(x),g(x)) are the indicial functions representing the change in the coefficient C a due to 

the unit step in a and q respectively. The indicial responses are functions of elapsed time, t - x, and 
are continuous single-valued functions of a(x) and q(x). The indicial functions approach steady-state 
values with increasing values of the argument t - x. To indicate this property, each indicial function 
can be represented as 


and 


c a a ( f - u o(t),$(t)) = C aa (oo; a(x),tf(x)) - ( t - x; a(x),q{x)) 


(2) 


c aq (t - x; a(x),q(x)) = C a? (oo; a(x),^(x)) - F Uq ( t - x; a(x),?(x)) 

where C a { oo; a(x),q{x)) is the rate of change of the coefficient C a with a(x) and q(x), evaluated at 
the instantaneous value of a(x) with q fixed at the instantaneous values of q(x). A similar definition 
applies for C a (oo; a(x),q{x)). The functions F a and F a are called deficiency functions. 

q a q 

When equations (2) are substituted in equation (1), the terms involving the steady-state parameters 
can be integrated and equation (1) becomes 
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C aW = C a( 00 ; a ( t U(0)-f F a a ( t-r l a(r),q( r)) a(r) dr 

o 

/ (3) 

- -J F a q ( f - F 4(t) dr 

0 

where C a (°°; a(t),g(t)) is the total aerodynamic coefficient that would correspond to steady flow with 
a and q fixed at the instantaneous values of a(t ) and q(t). 

Further simplification of equation (3) can be achieved by expanding the terms in this equation in 
Taylor series about q = 0, taking into account only linear terms and neglecting terms in q , qq and 
dq. Then 


Q(') = Q(°°; a ( t )^) +—C aq (<=o-, a(/),0) q(t) 


~f F a a ( l ~ T; a ( T )’°) a ( T ) dt 
0 


or in simple notation 


C a(0 = Q(«) + -^ C a (? («) < ?(0-J' F a a «(*)) a{t) dt 

0 

The deficiency function will be considered in the form 

F(t;a ) = h(t;a)a(a) 

where a(a ) is a polynomial in a, h{f, a) represents a sum of exponential functions 


(4a) 


(4b) 


(5) 


h(t;a)= e J 
j = 0 

and bj(a ) are again polynomials in a. For further analysis, however, only two forms of h(t;a) will be 
considered leading to the following deficiency functions 

F(t\a) = e~ blt a( a) (6) 
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and 


F(f,a) = e~ bl( ' a ' >t a{a ) (7) 

When the differential operator, D = dldt, is introduced and operator notation used, the convolution 
integrals with two forms of deficiency function can be expressed as 

J e~ b{ ( ?-T )a(ot(x))6t(x)uh: = -^—^-Da(t) (8) 

0 D + by 

and 


Je M a (* T ))( / T ^oi(x))d(x)uft = — — r Da(t) (9) 

0 D + b \\ a) 

Incorporating equation (8) into the operational form of equation (4) and recognizing that for one 
d.o.f. motion in pitch q = d results in 

C a (t) = C a (a) + ^C aq ( a)Da(t ) -^L Da (t) ( 10 ) 

Multiplication of both sides of equation (10) by ( D + b\) yields 


D C a (t) + by C a (t) = DC a (a) + byC a (a) + ~DC aq (a) Da(t) 


^- C a q (a)-o(a) 


Da(t ) 


( 11 ) 


Equation (II) can be considered as a postulated form for model identification, i.e. for model 
structure determination and parameter estimation. By examining equation (11), however, it is 
apparent that from measured time histories C a {t), a(t ) and their derivatives it is not possible to 

estimate explicitly parameter by and the remaining parameters in C a ( a), C Uq (a) , and a(a). To avoid 

this problem, it will be further assumed that 


a) C a (a) is known from static measurements, 
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b) C a ^(a) is estimated from small amplitude oscillatory data using techniques introduced in 
reference 8. 


Combining time histories C a (t ) with those of C a (a) and (f/V) C a (a) a(t ), a new variable y can be 
introduced as 

y( t ) = C a {t)-CaH-^C aq (a ) a(t ) (12) 

or, by using equations (4) and (6) as 


y(t) = -Je b, ^ r a(a(x)) d(x) dx 
0 

Equation (13) in operator notation will represent Model I as 


which is equivalent to 


y( '>— 7rf DaW 

D + b\ 


(13) 


(14a) 


y(t) + b x y{t) = -a( a) a (t) ( 14b) 

The second model considered incorporates dependency of the parameter by on the angle of attack. 
Equation (13) takes more general form (see Appendix A) as 


t -fbi (<*(!)) <€ 

y(t) = -Je T cf(a(T)) «("c) dx 

o 

Model II is then defined as 


or 


y(t) = 


a(a) 

D + 6j(a) 


Da(t) 


y(t ) + (?i(a) y(t) = -a( a) d (t) 


(15) 


(16a) 


(16b) 
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Motion in Roll and Yaw 


For a one d.o.f. motion in roll at a constant value of the angle of attack, o.q, the relations for the 
lateral aerodynamic coefficients are 


C a( t ) = C a(<P( t )>P( t . )); a = Y,i,OXn 

where the roll angle is related to the sideslip angle by the equation 

P = sin -1 (sin0sina o ) (17) 

The aerodynamic coefficients can be formulated as 

t 

C a {t) = C a (0)+r< c ap (t-r; P(r),p(r)) P(r) dr 
0 

r (18) 

+ C a p (t ~ V, P(t), p{ t )) p(r) dr 
0 

If the procedure illustrated for the motion in pitch is followed, equation (18) will be simplified as 

C a (t) = C a { cc; p(t),0) + ^C ap (co; p(t),o) p(t) 

' v . (19) 

-J F a p (f - T; p{r), o) P(t) dr 

o 

where C a (oo,(3(f),0) is the total aerodynamic coefficient that would correspond to steady flow at a 
fixed value of a and with (3 fixed at the instantaneous value of |3(r), and C a ^(°°,|3(f),0) is the rate of 
change of the coefficient C a with pit) evaluated at fixed value of a and instantaneous value of (3(0- 
F a is the deficiency function which might take the form of equation (5). 

Similarly, for a one d.o.f. motion in yaw at a constant value of a, the relations for the lateral 
coefficients are 


C a{ t ) = C a (ip(t),r(t))-, a = Y,£ or n 
where the yaw angle is related to the sideslip angle as 
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(20) 


(3 = sin ! (-sin^ cosag) 

The simplified model for the coefficients takes the form 

t 

C a (t) = C a {P) + ^C ar (l3)r(t)-J F a ^(t- r; fi(r)) ip(r) dr (21) 

o 

where the definitions of terms in (21) are similar to those for terms in equation (19). 

Model Identification 

Model structure determination and parameter estimation will be demonstrated on model equation 
(16) governing a one d.o.f. motion in pitch. Modifications to less complicated model (14) or models 
for a one d.o.f. motion in roll and yaw can be easily made. Substituting measured values at 
t[,i = 1,2 ,...,«, into equation (16b) gives 

f£(0 = -[^l («£:(()) JeW + «(«e(()) afO')] + £y(0 (22) 

where index E indicates the measured values, e y (i) is an equation error at time 0 -1) A t and At is the 
sampling interval. Equation (22) is the regression equation with the unknown parameters in 
polynomials b^a) and a(a). The mean values of these parameters can be estimated by a least squares 
technique. The parameter covariance matrix under the assumption of colored noise can be obtained 
from expressions in reference 11. In order to avoid differentiation of measured data an approach of 
reference 12 using modulating functions can be applied. 

For estimation of parameters in equation (22) a structure of both polynomials b\(a) and a(a) 
must be either known or determined from experimental data. The structure of a( a) can be selected 
from results of the small amplitude oscillatory data analysis as indicated in reference 8. If either the 
structure of b | (a) or structures of both terms, £q( a) and a{ a), are not known, a stepwise regression 
can be applied to model structure determination and parameter estimation (see e.g. ref. 13). 

The least squares parameter estimates can be updated by a maximum likelihood estimation 
method outlined in reference 11. The constraint equations are the state and measurement equations 
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of the form 

j(t) = -&l(a)y(/)-a(a)ci £ (r); y(t = 0) = y(0) (23) 

?£■(*■) = y{i) + v(i), i = l,2,...,n (24) 

where v(i) is the measrement noise at time (i - 1) At. 

In some cases the model for C a ^(a) may not be known. Then the parameters in b \ (a) and a(a) 
will be estimated for some a priori values of C a ^(a). Returning to equations (12) and (15), a new 
variable z(t) can be formed as 

t 

z(t) = C a (t)-C a (a)+j' e a(a) d(r) dr 

0 (25) 

= ^ C a q («)«(') 

When the measured values, and parameter estimates in b\ (a) and a(a) are substituted into (25) the 
regression equation is obtained as 

£ 

ZeO) = — d £ (/) + e,(/), i = l,2,...,n (26) 

Based on equation (26), the model structure of C a (a) can be determined and parameters in that 
model estimated. In the following step the parameters in b j(a) and a (a) can be estimated again, this 
time for the new model of C a ^ (a) and new values of y^d) computed from equation (12). This two- 

stage optimization procedure can be repeated until the minimum of the cost function for the 
maximum likelihood estimator is reached. A block diagram for the two-stage estimation procedure is 
presented in figure 1. 

Examples 

The procedure for identifying a nonlinear unsteady aerodynamic model of an aircraft subjected to 
one d.o.f. harmonic motion about one of its body axes is demonstrated in two examples. Both 
examples use data from pitch oscillations only. In the first example, the methodology is applied to 
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simulated data representing the pitching moment coefficient of a tailless aircraft. In the second 
example, wind tunnel data from a 10-percent-scale model of the F-16XL aircraft are used. A three- 
view of this model is shown in figure 2 together with some of the basic dimensions. Static and 
dynamic tests were conducted in the NASA Langley 12-Foot Low-Speed Wind Tunnel. A brief 
description of the test is given in reference 8. 


Example 1 

The purpose of this example is to demonstrate the feasibility of the algorithm to estimate 
parameters in the model with a given structure. In addition, the effect of measurement and modeling 
errors on the estimates will be investigated. The time histories of the pitching moment were computed 
from equations (4), (6) and (7), and data in table I for Model I and Model II. The expressions for 
by(a ) and a(a) were postulated as splines of the fonn 


a(a) 


= ciq + ct[a + fl2« 2 + a - a yj 2 
7=1 + 


(27) 


and 


4 / \ 

Ma) = b l0 + 2 fly (a - a y -) + (28) 

7=1 

where a j are knots and (a - aj) + are the plus functions defined as (a - ay) + = 0, for a j < a and 
(a - a j) + = a-a.j, for aj > a . 


The plots of the data in table I are presented in figure 3. Both the nominal value and amplitude of 
the angle of attack oscillations were selected as 35 deg, and the three frequencies of the oscillatory 
motion were 0.25, 0.50 and 1.00 Hz. The sampling interval was 0.01 sec. The variation of the 
pitching moment and its components with the angle of attack is shown in figure 4. The time histories 
of y and y are plotted in figure 5 for three cycles of each frequency. A zero-mean, Gaussian and 
white random sequence representing the measurement noise was added to the computed values of y. 
The variance of this sequence was defined by the signal-to-noise ratio (SNR). 

The effect of measurement noise on the ML estimates of b j and five parameters in a(a ) is shown 
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in table II. The increase of noise level, SNR changed from 40 to 20, resulted in expected increase of 
errors in estimated parameters, and in the fit error, y(v). The differences between parameter true and 
estimated values, however, remain within the 2o-confidence intervals. In table III the effect of 
modeling error in b i is demonstrated. Replacing the spline b^(a) by a constant led to a large fit error 
and large errors of parameters in a(a). The results in table III further indicate that even for correct 
structure of b\(a), the parameters were, in general, estimated with low accuracy. Finally, in table IV 
the results of two-stage optimization are shown. The parameter estimation started with an incorrect 
model for C m (a) by replacing the second-degree polynomial by a known constant. Then model 
structure determination and parameter estimation procedures were applied to identify a model for 
C m ^(a) in regression equation (26). As indicated in table IV, after three iterations the identified 

model for C m ^(a) was very close to the true one. The remaining parameter estimates were also close 
to their true values. 

Example 2 

The measured static and oscillatory data used in this example are shown in figure 6 as C L ( a), 
C m ( a), C L ( a; a 0 , a A , k), and C m (a; a 0 , a A , k ) where a 0 = 35 deg, a A = 35 deg and k = 0.034, 0.057, 

0.1013 and 0.1350. For the wind tunnel speed V = 17.52 m/sec and c = 0.753 m, the corresponding 
frequencies were /= 0.25, 0.42, 0.75 and 1.00 Hz. Each of the four time histories of the oscillatory 
data were comprised of three cycles with the sampling rate of 100 Hz. The time histories of 
measured data were obtained as the average values from five repeated runs at the same amplitude and 
frequency. The variability of averaged data in cycles was, in general, very low. Some scatter 
appeared in the stall region of the pitching-moment coefficient as can be seen in figure 7, where the 
data from three repeated cycles are shown. The analytical forms of static data were obtained by 
fitting the measured C L (a ) and C m (a) curves. For the a priori values of two damping terms, C L (°°;a) 

and C m (°°;a), the estimates from small-amplitude oscillatory data of reference 8 were used and 
reformulated as 
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and 


for a <20 deg 
for a < 20 deg 


C Lq ( co;a) = -0.424 

= -2.0 + 5.7a - 3.4a 2 


C m (oo;a) = -1.245 - 0.3806a + 1.5557a 2 

The models for the polynomials a(a) and b ( (a) were postulated as polynomial splines given by 
equation (27) and (28) with two knots in each expression. The variable y E was computed from 
equation (12), its derivative was obtained by numerical differentiation. 

After two iterations of the two-stage optimization algorithm the identified models for the 
polynomials a(a ) and b \ (a) were 



6[(a) = b 0 + b^a + (30) 

7 = 1 


for the coefficient C L (a ) and 


c?(a) 


= 



for the coefficient C m (a). 


6i(a) = Z?io 


(31) 

(32) 


The ML estimates of model parameters in equations (29) to (32) and their standard errors (Cramer- 
Rao bounds) are summarized in table V. The standard error of estimated parameters varied between 
4 to 11 percent indicating possible identification problems for some parameters in the model. The 
plots of polynomials a (a) and 6|(a) are presented in figures 8 and 9. The a priori and estimated 
values of parameters in polynomials representing the variation of the damping terms C rn ^ with the 
angle of attack are shown in table VI. The identified model had the same structure as its a priori 
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counterpart, and the accuracy of the estimated parameters was between 3 to 13 percent. As pointed 
out in table VI the a priori model for (<»;a) was not updated because of a small contribution of 

the Ci term to the lift. The identified final models fit the measured data very well at all frequencies. 

An example of measured and estimated coefficients is given in figure 10 for the reduced frequency 
k = 0.057. 

The identified models were also assessed by their prediction capabilities. The predicted time 
histories of C L and C m were computed from equation (4b) for selected amplitude and frequency of 
the oscillatory motion or for the ramp input in the angle of attack at different rates. A comparison of 
measured and predicted coefficients C L {a) and C m ( a) for two different amplitudes and similar 
frequencies is given in figures 11 and 12. Figures 13 and 14 present a comparison of the same 
coefficients for two different ramp inputs versus a. The same data in the form of time histories are 
shown in figures 15 and 16. The results in figures 11 to 16 indicate that the identified models are 
good predictors for the lift coefficient, while some discrepancies between measured and predicted 
data can be seen in the pitching-moment oscillatory data with the amplitude of 20 deg and the ramp 
data. 

Concluding Remarks 

Aerodynamic equations with nonlinear unsteady effects were formulated for an aircraft in a one- 
degree-of-freedom large amplitude motion about each of its body axes. The corresponding 
aerodynamic models were expressed in the from of indicial functions. The model formulation 
separated the resulting aerodynamic forces and moments into static terms, purely-rotary terms and 
unsteady terms. The unsteady term in the model for a pitching motion was modeled as a product of 
an exponential function and a polynomial in the angle of attack. For model identification from 
experimental data it was assumed that the static and purely-rotary terms were known. The model 
identification procedure developed combines stepwise regression and maximum likelihood 
estimation. In cases when the a priori information about the rotary term is in doubt, a two-stage 
optimization algorithm which can identify both the unsteady and rotary terms were proposed. 
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The identification scheme was applied to wind tunnel oscillatory data in pitch in two examples. 
The first example used the simulated data for a tailless aircraft and the second used wind tunnel 
oscillatory data from the F-16XL aircraft. The results from both examples indicated that 


1. the two-stage optimization algorithm can converge to maximum likelihood estimates; 

2. the accuracy of estimated parameters can be severely degraded by modeling errors; 

3. the identified model from experimental data fit the data well, however, the accuracy of some of the 
estimated parameters was rather low, around 10%; 

4. the identified model was a good predictor for oscillatory data and data with ramp input. 
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Appendix A 

Integral form of Model II 

The following differential equation is considered 

y(t) + bi(a) y(t) = - a(a) d(t) (Al) 

where b^a) and a(a) are polynomials in a. After multiplying each side (Al) by the exponential 
term 


fbi{a 


and rearranging, the following relationships is obtain 


d ( 

— y(t)e° 

at 


= - a(a(t )) a (t)e 




\ ! 
Integration of both sides of (A2) results in 



' fb ,(«(§)) 

-J' a(a(r)) d(r)e° dr 


or 


y(t) = -fe 0 ' £7(a(-c)) oi(t) dr 

o 

Equation (A4) leads to the final form expressed as 


t -vMg))dg 

y(t) = -Je T a(r) dx 

o 


(A2) 


(A3) 


(A4) 


(A5) 
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Table I. Characteristics of a tailless aircraft used in generating oscillatory data 


C m (a ) = 0.42a + 0.34a 2 - a 3 + 0.40a 4 
C mq (a) = -2 + 0.8a 2 

a(a) = +0.28a - 3.2a 2 + 8(a - 0.4363) 2 
-7.4(a - 0.9599) 2 

Model I: 
b 1 = 2.5 
Model II: 

b x {a) = 2.5 - 5.73(a - 0.349) + + 5.73(a - 0.5236) + 
+5.73(a - 0.827) + - 5.73(a - 1 ,0472) + 

— = 0.02131 sec 
1 / 


Table II. Effect of measurement noise on estimated parameters. 
Simulated data. Model I. 


Parameter 

True Value 

Estimate 

SNR = 40 

SNR = 20 

b \ 

2.5 

MiBlS 


a \ 

.28 

.28 

(.010) 

mm 


-3.2 

-3.21 

(.022) 

-3.17 

(.046) 


8.0 

8.0 

(0.44) 

7.91 

(.092) 

A 2 

-7.4 

-7.6 

(.20) 

-6.8 

(.42) 

s(v) 

— 

.0041 

.0085 


Note: numbers in parentheses are Cramer-Rao bounds on standard 
errors 
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Table III. Effect of modeling errors in b | (a) on estimated 
parameters. Simulated data, Model II, SNR = 40. 


Parameter 

True Value 

Estimate 

Model II 

Model I 

^10 

2.5 




-5.73 

-5.9 

U2) 

— 

b 2 

5.73 


— 


5.73 


— 

B 4 

-5.73 

-5.9 

(.39) 

— 

a l 

.28 


SmmMSwm 

“2 

-3.2 

-3.16 

(.029) 

-4.25 

(-047) 

Ai 

8.0 



a 2 

-7.4 

-7.0 

(.25) 


s(v ) 

— 

.0045 

.0087 


Note: numbers in parentheses are Cramer-Rao bounds on standard 
errors 
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Table IV. Effect of modeling error in C m (a) on estimated 

<? 

parameters. Simulated data. Model I, SNR = 40. 




Estimate 

Parameter 

Tme 

Value 


C m q (a) 
estimated 

b \ 

2.5 

■■ 

2.464 

(.0033) 

a \ 

.28 

.55 

(.022) 

.32 

(.010) 


-3.2 

-3.62 

(.050) 

-3.26 

(.023) 

^1 

8.0 

8.49 

(.099) 

8.08 

(.045) 

A 2 

-7.4 

-8.1 

(.45) 

« 

s(v) 

— 

.0092 

.0042 


Note: a) number in parentheses are Cramer-Rao bounds on 
standard errors 


b) initial value C m =-1.4 

1 

2 

estimate after four iterations: C,„ = -1.939 + 0.804a 

(.0043)(.0072) 

true model: C„, = -2 + 0.8a 2 
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Table V. Maximum likelihood estimates of parameters ina(a) and 
^(a) polynomial splines. 


Parameter 

C L 

C L knot 
location 
(degrees) 

c m 

C m knot 
location 
(degrees) 

% 

— 

— 


— 

a \ 

9.8 

(.56) 

— 

5.1 

(.43) 

— 

°2 

-31.0 

(1.4) 

— 

-9.7 

(.69) 

— 

A l 

50.0 

(2.3) 

20 

18.0 

(1.2) 

26 

A 2 

-22.0 

(2.5) 

47.5 

-11.0 

(1.3) 

46 

b 0 

12.6 

(.45) 

— 

7.1 

(.36) 

— 

b l 

-16.1 

(.78) 

— 

— 

— 

B i 


45 

— 

— 

*2 

-32.0 

(4.5) 

55 

— 

— 

con 

-0.14 

(.018) 

— 

0.07 

(.010) 

— 

5(V) 

.026 

— 

.012 

— 


Note: (a) numbers in parentheses are Cramer-Rao bounds on 
standard errors. 

(b) con is a constant added to the state equation. 
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Table VI. A priori values and least-squares estimates of parameters in C Lq {^;a ) 


Parameter 

with 

Q^ 00 ; «) 

C mq(^ «) 

a < 20 deg 

a > 20 deg 

a priori 

estimated 

a 0 

-0.424 

-2.0 

-1.245 

■a 

a 

0.0 

5.7 

-0.381 

-0.76 

(.098) 

a 2 

0.0 

-3.4 

1.556 

n 


Notes: (a) numbers in parentheses are Cramer-Rao bounds on standard errors. 

(b) C mq parameter estimates are obtained after two iterations. 

(c) C [q was not updated due to its small contribution to C L . 
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two stage optimization 


fyes 


yE 



Eqn. A: z(t) = [fy C a q tfW 

Eqn. B: y(t) = C a {t)-Cj~\a,Qi)-^C a (~;a, 0) <?(f) 

Figure 1. Block diagram of model identification using stepwise regression (SR) and a maximum likelihood (ML) 
estimation. 
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Figure 2. Three- view sketch of F-16XL model. 
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a, deg 


Figure 4. Pitching-moment coefficient and its components in steady oscillatory motion at /= 0.5 Hz. Simulated 
data, (a) Model I, (b) Model II. 
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